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On  Dynamic  Plastic  Mode  Form  Solutions 


by  ?.  S.  Symonds^  and  C.  T.  Chon^ 

Abstract 

The  paper  discusses  the  sippificance  and  calculation  of  dynamic  plastic 
mode  form  solutions  for  small  deflections  of  structures  of  ripid-nerfectly 
plastic  materials  subjected  to  load  systems  of  fixed  distribution  and  magnitude. 
These  solutions  have  separated  form,  with  velocity  the  product  of  a scalar 
function  of  time  by  a vector-valued  function  of  space  variables.  The  relation 
is  shown  on  the  one  hand  to  mode  form  solutions  for  a structure  of  viscoplastic 
material,  and  on  the  other  hand  to  limit  load  solutions  for  those  of  perfectly 
plastic  behavior.  Numerical  examples  are  piven  for  circular  plates  of  material 
obeying  Tresca  and  Mises'laws. 


1.  Introduction.  "Mode  form"  solutions  of  the  equations  poveminp  the 


■ 


nlastic  response  of  a structure  to  d3mamic  loads  are  "particular  ince^als"  in 
separatee  foirni:  the  velocity,  for  example,  is  the  product  of  a function  of 

space  variables  by  a function  of  time.  In  this  paper  we  are  concerned  with 
solutions  of  this  type  which  hold  throuf^hout  the  motion,  which  we  call  "permanent" 
mode  form  solutions.  They  satisfy  at  all  times  the  stated  field  equations  of 
dynamics,  kinematics,  and  material  behavior,  together  with  loading  and  boundary 
fixing  conditions  [1].  As  will  be  seen,  such  solutions  exist  for  certain  ideal- 
izations of  material  behavior  and  loading,  and  when  geometry  changes  are  neg- 
lected in  the  governing  equations.  They  are  to  be  distinguished  from  "instantan- 
eous" mode  form  solutions  which  satisfy  the  field  equations  only  in  a special  sense 
[2  - 6]:  regarding  the  displacement  field  as  instantaneously  fixed,  and  a scalar 

magnitude  of  the  velocity  field  also  as  fixed,  an  instantaneous  mode  fonn  solu- 
tion has  acceleration  and  velocity  fields  over  the  structure  that  are  identical 
apart  from  a scalar  factor.  The  common  shape  function  satisfies  the  field 
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equations  appropriate  to  the  stated  i nst ant  cine ous  displacement  function  and 
velocity  mappitude. 

Both  tvnes  of  mode  form  solutions  have  been  used  to  obtain  relatively 
simple  approximate  solutions  of  problems  of  structures  loaded  impulsively, 
i.e.,  by  very  short  hiKh-intensity  pressure  pulses  whose  effect  can  be  spec- 
ified by  a field  of  initial  velocities,  with  initial  displacements  taken  as 
zero.  Structures  both  of  conventional  plastic  (time  independent)  and  of  rate 
dependent  plastic  (viscoplastic)  material  behavior  have  been  treated.  Perm- 
anent mode  form  solutions,  when  they  exist,  provide  extremely  simple  solutions 
of  initial  motion  problems  that  are  appropriate  for  conditions  in  which  geo- 
metry changes  are  negligible.  However,  in  certain  technically  important 
cases  the  finite  deflections  have  a large  influence  on  the  response  and  must 
be  included.  Approximate  solutions  for  such  cases  have  been  constructed 
from  a sequence  of  instantaneous  mode  form  solutions,  appropriate  to  successive 
stages  in  the  response  [5,  6].  In  both  types  of  application  the  specified 
initial  velocity  conditions  are  not  satisfied  exactly,  but  only  in  a "least 
squares"  sense.  This  starting  point  also  provides  a measure  of  error,  useful 
at  least  for  comparing  two  such  solutions. 

Fundamental  properties  of  dynamic-plastic  structural  motion  underlie 
these  methods  [1-4].  The  response  of  a structure  to  dynamic  loading,  how- 
ever started,  tends  toward  a mode  form  solution.  Both  types  of  mode  solu- 
tions render  a certain  functional  an  extremum.  Linear  equations  of  dynamics 
and  kinematics  are  presumed  for  these  properties,  which  hold  for  a wide  class 
of  plastic/viscoplastic  behavior. 

Some  of  these  properties  are  illustrated  by  a very  simple  structural 
model  with  two  degrees  of  freedom  [5,  7].  As  shown  in  Fig.  1(a),  the  model 
simulates  a simply  supported  beam  by  two  mass  particles  connected  by  rigid 
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massless  rods,  with  loints  at  the  two  masses.  The  loads  are  forces 
applied  to  the  masses.  Deformation  can  occur  only  at  the  ioints,  and  the 
material  behavior  is  specified  in  terms  of  moment- -rotation  rate, etc.  character- 
istics. For  riEcid-perfectly  plastic  behavior,  the  characteristic  curve  is 

shown  in  Pip,  1(c);  fully  plastic  moments  ± M are  associated  with  rotation 

o 

rates  of  correspondinp;  sign  and  arbitrary  magnitude.  The  rotation  rates 

• • * • 

'*’1’  "^2  ®re  related  to  the  transverse  velocities  u^,  u^  bv 


(la,b) 


The  equations  of  motion  of  the  two  masses  can  be  written 


Pi  - Yi  = Gu^ 


P2  - V2  = "-2 


(2a, b) 


where  G is  the  mass  of  each  particle  and  Y^,  Y^  are  the  structural  re- 
action forces;  for  small  deflections,  with  notation  as  indicated  in  Fig.  1(d), 


= . L.  (-K^  . 3«^) 


(3a, b) 


We  take  the  bending  moments  as  satisfying 


-M  < M < M (4) 

o - a - o 

vrfiere  a = 1,  2.  Considering  these  inequalities,  the  locus  of  forces  Y^ 

such  that  either  or  is  equatl  to  ± is  shown  as  curve  ABCD  of 

Fig.  2(a).  By  definition  this  is  the  "yield  surface"  of  the  structure.  For 

P = P^  = Y^  , unlimited  slow  deformation  can  occur,  while  load  states  P 
a o a a 

inside  the  curve  can  be  supported  without  deformation;  none  occurs  unless  due 
to  previous  loading  history.  Kinematic  relations  obey  the  standard  normality 
rule  of  yield  surfaces;  here  the  "strain  rate"  vectors  of  the  structure  are  the 
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velocity  vectors  = (u^,  ) . For  example,  side  AB  of  the  structure  yield 

curve  corresponds  to  M,  = M , -M  < M.  < M . For  any  load  tx)int  not  at  either 
^ 1 0 0-2-0 

of  the  corner  points  A,  B the  velocity  diagram  has  the  shape  shown  in  Fif». 

2(b)  with  Uj^  = Bu^  , since  = 0 for  jM^I  • The  velocity  vector, 

as  indicated  in  Fig;.  2(a),  is  pemendicular  to  the  line  AB  , whose  eouation  is 

3Y,  + = 4M  /£. 

12  o 

The  discussion  above  assumes  static  behavior.  Now  we  consider  dynamic 
problems,  the  structure  havinf^  non-zero  accelerations  u^  as  well  as  veloci- 
ties u^  . The  dynamical  equations  (2)  mav  be  written  in  vector  form  with 
the  dimensionless  variables  as 


Y = u 


(5) 


where 


i:  = 


V = ir(Vj.Vj) 

o 


Gl  /•  " . 

y ( u , Uj  ) 


(6) 


A typical  dynamic  case  is  shown  in  Fig.  2(a)  for  a load  vector  P'  which  can- 
not be  supported  statically.  The  instantaneous  stress  vector  Y and  velocity 
vector  u must  be  related  through  the  normality  rule.  Suppose  the  velocity 

A 

vector  has  the  direction  of  vector  u'  . The  stress  state  vector  must  cor- 
respond to  the  comer  point  B for  this  direction,  which  lies  between  the  nor- 
mals to  the  adjacent  sides  at  point  B.  The  acceleration  vector  u is  required 

by  Eq.  (5)  to  have  direction  and  magnitude  as  shown.  Evidently  this  is  not  a 

**  • 

mode  form  motion,  since  the  vectors  u and  u are  not  parallel.  However,  it 
can  be  seen  that  if  the  load  P'  remains  constant  the  velocity  vector  must 
change  with  time  until  its  direction  and  that  of  the  acceleration  vector  become 
normal  to  side  AB  of  the  yield  diagraf’'.  This  "mode  form"  motion  is  indicated 
by  the  dashed  vectors.  This  example  illustrates  the  convergence  pr>opertv  on 
v/hich  the  mode  approximation  technique  is  based  [1]. 
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Several  different  types  of  mode  form  solution  are  illustrated  in  Pip;. 

2(b).  They  have  the  common  feature  that 

u = A(t)u  (7) 

a a 

where  A is  a scalar  factor.  We  consider  forces  in  fixed  ratio 

to  each  other.  Then  the  vector  P can  be  written  as 

P = XP  (8) 

where  A is  a scalar  load  factor  and  P is  a constant  nominal  load  vector. 

In  the  examnle , P = (1,  - 0.5).  Fif^ure  2(b)  shows  a number  of  cases  cor- 
responding to  nosit ive  values  of  A . At  the  plastic  collapse  or  limit  load 
with  A = Aj^,  the  acceleration  is  zero  but  the  velocity  is  nonzero;  this  is  case 

*( 3 ) *( 3 ) 

3,  with  Y = P . The  other  extreme  is  that  of  P = 0,  shown  as  case  1 with 
= -u^^\  Here  it  is  implied  that  some  previous  impulsive  (or  other)  load- 
ing history  has  resulted  in  the  mode  pattern  Fig.  2(b)  with  u^  = Su^  , so  that 
the  velocity  vector  has  the  direction  of  the  normal  to  side  AB  and  is  par- 
allel to  -u  . Note  that  the  initial  conditions  and  nrior  loading  history  are 
relevant  only  in  the  sense  of  causing  a particular  mode  form  to  occur  rather 
than  another  possible  one.  Evidently  for  A < Aj^  the  mode  solution  is  not 
unique;  for  A = 0 there  are  eight  modes,  with  stress  vector  either  perpendi- 
cular to  one  of  the  four  sides  or  at  one  of  the  four  comer  points  A,  B,  C,  D. 

(2 ) 

The  case  marked  as  case  2 has  A ^ mode  is  indicated  by  the 

velocity  vector  u'  . In  this  case  there  are  two  possible  mode  solutions, 
the  alternative  mode  shape  having  velocities  proportional  to  ( -1,  -3),  with 
acceleration  vector  perpendicular  to  side  BC  of  the  yield  curve,  as  indicated 
by  dashed  vectors  in  Fig.  2(b). 


f 
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The  case  marked  as  case  4 has  X > X^  , so  that  the  load  state  point 

lies  outside  the  yield  curve.  For  the  example  shown  there  is  only  one  mode 

r.(4)  . . . 

solution,  the  acceleration  u having  the  same  direction  as  the  velocity 

• ( 4 ) , , 

u . However,  note  that  for  X = X„  the  stress  state  ooint  is  at  comer  B. 

- D 

Case  5 has  X^^^  > \ , The  velocity  vector  can  have  anv  direction  between 

B 

the  normals  to  the  adjacent  elements  at  B,  so  that  the  common  direction  of 
the  acceleration  and  velocity  vectors  changes  as  X increases.  For  X -+■  * 
both  anproach  the  direction  of  the  load  vector  P . 

2.  Modes  in  a Circular  Plate.  To  illustrate  the  reouirements  for  mode  form 
solutions  of  a continuous  structure  whose  material  exhibits  either  rigid- 
perfectly  plastic  behavior,  or  rigid-viscoplastic  (strain  rate  dependent)  be- 
havior, we  consider  a circular  plate  of  radius  R , with  a radially  symmetric 
pressure  p(r,t)  applied  transversely.  As  above,  we  consider  only  the  plate 
in  its  initial  configuration,  i.e.,  we  are  not  considering  Questions  of  finite 
deflections.  The  equation  of  motion  is 

If  [lF^^”r^  - Mg']  = - pr  + prw  (9) 

Fixing  conditions  are,  for  example 

w(R,t)  = w’(R,t)  = 0 (10) 

Initial  conditions  may  also  be  stated  as 

w(r,0)  = 0 ; w(r,0)  = w°(r)  (11) 

In  these  equations,  r and  0 are  polar  coordinates  in  the  middle  surface 

and  t is  time;  w(r,t)  is  the  transverse  displacement,  w'  = 3w/3r  the 

• "22 
Slone,  w = 3w/3t  the  velocity,  and  w(r,t)  = 3 w/3t  the  acceleration; 

M^  , Mg  are  radial  and  circumferential  bending  moment  per  unit  length,  res- 
pectively; and  p is  the  mass  per  unit  area  of  middle  surface. 
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The  strain  rate  quantities  relevant  to  the  plate  in  bending  are  curva- 
ture rates  ic  , ic.  defined  so  that  the  work -dissipation  rate  equation  is 
r 0 


rR 


(p  - pw)  w rdr 


(M  IC  + M.ic.)  rdr 
r r 0 0 


(12) 


Eauations  (9)  and  (12)  are  consistent  with  the  following  ^orms  of  the  curvature 
rates  in  terms  of  the  plate  velocity 

1 


IC  = - w’ 

r 


ic„  = w' 

0 r 


(13a ,b) 


We  will  write  constitutive  equations  suitable  either  for  rigid-perfecrtly 
plastic  material  or  for  one  with  strain  rate  sensitivity.  Both  behaviors  can 
be  treated  by  writing  constitutive  equations  that  are  generalizations  of  the 
following  form  for  simple  tension: 


sgn  o 


(14) 


This  homogeneous  viscous  form  can  represent  adequately,  over  a range  of 


strain  rate,  the  more  realistic  inhomogeneous  form 

in 


1^1  -1 


sgn  c for  jo  I > o 


(15a) 


e = 0 for  I o|  ^ Og  (15b) 

This  latter  form  is  capable  of  close  representation  of  curves  for  stress  as 
function  of  strain  rate  at  constant  strain  [8];  it  is  assumed  that  strain  rate 
history  effects  are  absent.  The  viscous  form  of  Eq.  (14)  involves  no  yield 
condition.  We  need  to  write  stress-strain  rate  relations  which  generalize 

Eq.  (14)  to  appropriate  forms  for  a thin  plate  in  bending,  in  terms  of  bending 

• • 

moments  M M„  and  curvature  rates  k_  » • 

■p  U r (7 


We  write  for  brevity 


8 


in 


' = M /M’  , C = ic  /<  where  a = r or  9,  /**»  *'o  ” 


a a o ’ a a o 
Then  *we  can  write 


c - r"'  II.  - ^ 

^ a™'  " am’ 


3m'  9m’ 

a a 


a6a  ) 


where 


F = / m'^-  m'm'  + 

r r*  u “ 


(16b) 


ir +1 


(l6c) 


and 


, . -1/n’  3G  _ 3« 

m - o - r 


(l7a) 


where 


G = ? /il^  i 

/r/ 


(l7b) 


n'  _l+l/n' 
n'+l 


(17c) 


These  forms  are  derived  by  assuming  a sandwich  plate  model  of  thickness 
h = H/2,  H being  the  thickness  of  the  uniform  plate.  Each  sheet  of  the 

sandwich  plate  is  in  plane  stress,  and  has  principal  strain  rate  components 

• • 

t , c-  given  in  terms  of  stresses  o , a.  by 
r o r o 


a 

a 


O’ 

o 


(18a, b) 


V * r~^  5"  2 /5  T":  ^ 

where  f = / % - ^ 

and  a = r,  9 . The  forms  of  Eqs.  (16 » 17)  are  capable  of  close  representation 
of  strain  rate  behavior  over  appropriate  strain  rate  ranges.  If  n*  is  made 
very  large  while  o^  is  taken  as  the  static  yield  stress  o^  , we  obtain  rigid 
perfectly  plastic  behavior. 
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We  look  for  a mode  form  solution  of  Eq.  (9)  satisfying  the  edge  conditions 
(10),  kinematic  relations  (13),  and  constitutive  equations  (15),  (17).  We 
write  the  velocity  as 

w = w^(t)*(r)  (19  ) 

and  take  i^(0)  = 1,  so  that  w^(t)  denotes  the  midpoint  velocity.  The  cur- 
vature rates  are  then 


K 

V 


w*k„ 
* 0 


= w*(-  ^V) 


(20) 


The  bending  moments  can  also  be  written  in  separated  variable  form,  because  of 
the  homogeneity  of  u)(E^),  as 


r 0 0 * 3k 


u w , * 1/n ' 3(1) 

M„  = K M'w* 

9 o o * 3k 


(21a,b) 


e 


where  )jj(r)  = d)[k  (r)].  The  equation  of  motion  then  can  be  written  as 
a 


3 , 3(2  , 3u) 

^ 3k  ^ 

• -1/n' 

= 

" Kt)p(r) 

r 0 

— — 

(22) 

Here  we  have  further  assumed  that  the  load  pressure  has  a constant  distribution 
pattern  over  the  plate,  and  have  written 

p(r,t)  = A(t)p(r)  (23) 


For  a rate  sensitive  material  Eq.(22)  can  be  "separated"  only  if  the 
pressure  distribution  p(r)  is  the  same  as  the  mode  shape  function  ♦(r) 
to  within  a constant  factor.  We  conclude  that  apart  from  this  highly  special 
case,  mode  form  solutions  cannot  exist  when  the  plate  is  subjected  to  a pres- 
sure loading,  even  for  the  homogeneous  viscous  behavior.  If  p(r)  = 0,  separation 
of  variables  is  possible  since  one  can  write 
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1 

pr* 


9 I 9 , 9(1)  > 9(1) 

9F  ^ ^ " 9E 

'—  r 0 


1 


K M' 
o o 


w.-. 

1/n ' 

w.-. 


A 


(24) 


wherv)  A is  a constant.  It  is  seen  that  in  this  case  the  midnoint  velocity 

• • • * o 

decreases  from  its  initial  value  w.,.  according  to  the  exoression 


w.^(t)  = 


1- 


c Cl  1 1 * » I A*  n -1 

( w... ) - ( 1 - — , ) ic  M ' At 

n ' o o 


(25) 


The  share  function  'ti(r)  must  be  obtained  so  as  to  satisfy  the  ordinary  dif- 
ferential equation  obtained  by  setting  the  left-hand  side  of  Eq.  (21)  equal  to 
-A.  It  may  be  solved  by  iterative  methods  [6,  8]  without  difficulty.  When 
$(r)  has  been  obtained,  the  magnitude  A can  be  calculated  from  the  energy- 
dissipation  rate  Eq.(12),  which  takes  the  form 


• ** 

fR 

- 2 • . 

p4i  rdr  = x M' 
0 o 

• 

0 

where  = 

=/r 

)2  t 

r 

C X(  ((i ) ] iP  rdr 
o 


(26) 


because  the  dissipation  rate  can  be  expressed  in  terms  of  the  homogeneous  form 
(D  .using  Euler’s  formula. 

For  rigid-perfectly  plastic  behavior  we  put  ~ n ' . Then  sol- 

utions in  separated-variable  form  exist  not  only  for  the  impulsive  loading  case 
with  p(r)  = 0 , but  for  arbitrary  p(r)  , provided  X(t)  = constant.  For  this 
material  the  stress  vector  depends  only  on  the  direction  of  the  strain  rate 
vector,  not  on  its  magnitude,  so  that  in  a mode  form  solution  the  stresses  are 
constant.  Thus  if  the  load  factor  X(t)  = X = constant,  the  acceleration  is 
constant  and  given  by 


Mm 
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fR 


'TP. 

jpiji  rdr 


p6rdr 


(M  k + M„k„)rdr 

TV  a 7 


f27) 


The  acceleration  evidently  can  be  ne, {native  (as  in  the  impulsive  loadinp;  case 
where  X = 0) , or  positive;  zero  acceleration  is  the  case  of  "plastic  collapse" 
in  conventional  plastic  limit  analysis  theory.  This  ranp;e  of  possibilities 
has  already  been  illustrated  by  the  two-mass  model. 

3.  Numerical  Results  for  a Circular  Plate.  We  consider  henceforth  a plate  of 
ripid-perfectlv  plastic  material  subjected  to  a uniform  pressure  p of  various 
macpitudes.  The  plate  has  uniform  mass  density  p per  unit  .\rea,  and  constant 
thickness  H.  IntegratinR  and  assuming  the  mode  form  En.  (19),  the  equation 
of  motion  Eq.  (9)  can  be  written  in  nondimens ional  form,  with  m = M /K  , 
m.  = Mq/M  , and  putting  r/R  r,  0 < r ^ 1,  as 

7 7 0 


) 


‘ 7 K ' 


M 

I o 


(fcrdr 


(28  ) 


The  work-energv  rate  Eq.  (24) becomes 


rl 

(^rdr  - 
- o 


PW^jR2 

M 

o 


rdr  =-  i (m4)"+m 
j I' 

‘ 0 


0 


(29) 


At  a "hinge  circle"  the  slope  nay  be  discontinuous,  and  the  right  hand  side 
must  then  include  a term 


-r.m  (r, ) At'(r  ) 
i r 1 i 

1 ■ 


- r.m  (r,  )[(ti'(r  + 0)  - 4>'(r  - 0)]  (30) 

1 r i 1 i 


for  a hinge  circle  at  r = r 


We  have  boundary  conditions  at  the  edge  r = 1 


1? 


M 


Simplv  supported  edpe  = <t>(l)  = 0,  m^(l)  = 0 (31  a) 

Clairoed  edpe  =(t(l)  = 0,iti'(l)  = 0 (31b) 

At  the  center  of  the  plate  r = 0,  m (0)  = m (0)  bv  S’/rtimetr"/.  V.’e  adopt  the 
normalizing  condition  <>(0)  = 1,  which  merely  defines  w*  as  the  midpoint 
velocity  of  the  plate. 

(A)  If  we  adopt  Tresca's  yield  condition  and  the  associated  flow  rule 
(Fiy.  3),  simple  closed  form  solutions  can  easily  be  derived.  The  analysis  is 
an  obvious  extension  of  that  for  the  limit  analysis  problem,  which  corresponds 
here  to  putting  w^^  = 0.  For  the  simply  supported  edge  case,  the  stress  states 
lie  on  segment  AB  of  the  yield  curve,  and  all  equations  can  be  satisfied  with 
velocity  profile  either  conical  in  shape  as  in  Fig.  4a,  or  with  shape  of  a 
truncated  cone  as  in  Fig.  4b.  The  conical  profile  holds  for 

0 ; ^ < 2 <32) 

o 

and  the  relation  between  load  and  acceleration  is 


2 2 
pw^r  . pR 

12M  ■ 6M 


( 33) 


o o 

When  p is  larger  than  the  upper  limit  set  by  Eq.  (32),  a violation  of  the 
yield  condition  occurs  at  r = 0 ; thus  the  truncated  cone  field  of  Fig.  4b  is 


required.  This  holds  with 


o o 

while  Eq.  (34)  gives  the  parameter  n defining  the  (constant)  shape  of  the 


response  field  for  large  pressures  (more  than  twice  the  limit  load  magnitude). 


For  a olatnped  plate  the  stress  states  lie  on  seprnent  AB  of  the  yield 


curve  for  0 < r < C , but  lie  on  BC  for  f;  < r < 1,  in  order  to  satisfy  the 

edve  conditions  Eos.  (31b).  The  strain  rate  vector  (r  , r.)  = Wj^  ( 

r r ••  r 

TTiust  be  oarallel  to  the  -m  axis  at  r = 1,  since  4'(1)  = 0 , hence,  the 

r 

flow  rule  requires  the  stress  state  point  at  r = 1 to  be  point  C of  the 
Tresca  diasrram.  Thus,  if  the  pressure  is  not  too  lar^e,  the  resnonse  field 
has  a conical  shape  for  r < (I  and  a loj^arithmic  one  for  E < r < 1,  as 
sketched  in  Tip.  4C.  The  mode  function  is  ptiven  by 

0<r<C:4>  = l+  .-^Tln^  g-l)  <35a  ) 

C < r < 1 : t (35b) 

The  pressure  and  acceleration  are  expressed  in  terms  of  K as  follows: 


pP^  _ 14^^in  c - 4f;^(in  r.y  - * G 

’ 2 ? ' *7 

o ^ (2r/ in  ♦ f In  ^ ♦ 3> 


o w . R 2 2 9 

° * _ 14r/in  r - 4r/(in  r)‘  _ lor/  - 6 In  F.  + 6 

b ' ~2  ? ' 5 

o r (?r  In  f.  + 6 InF,  - 3f,  + 3) 


provided  0.6201  < C < 0.8055 
The  important  special  cases  are 


Free  deceleration,  p = 0 : C = 0.6201  ; 


Limit  load  pressure,  w.  = 0 : p = 0.7300  ; 


dw^R 


PlR2 


= - 23.58 


= 11.26 


( 36a ) 

(36b) 

(36c) 

(37a) 

(37b) 


A violation  of  the  yield  condition  at  r = 1 will  occur  for  pR  22  49 
corresponding  to  ^ > 0.8055,  if  the  above  mode  field  is  retained.  The  mode 
shape  for  large  pressure  magnitudes  has  constant  velocity  in  an  interior  region 
as  sketched  in  Fig.  4d.  The  mode  function  is 


14 


0 ( r ( n : 

r)  f V f E,  : 


* = 1 

. _ V + E In  E - E 
^~E  In  E-E  + n 


(38a) 

(38b) 


e $ r ^ 1 ; 


C In  r 

^ In  C - 5 + 0 


(38c) 


The  load  nressure  and  acceleration  are  related  to  n and  E by 


9 — r.2 

pR  _ _ 12C  (C  In  C - C + n) 

M ■ M " 4 3 ^,3  ^4 

o on-  2En  * 2E  n - E 


(39a) 


(1  - In  On"^  - 2^n^  + (5?^  - 3^  - 2E^  In  On  - (E^  + 3^^  In  C)  = 0 (39b) 


nrovided  n i 0,  ^ ^ 0.8055,  5 22.49 

O 


(39c) 


These  results  for  the  static  limit  pressure  and  for  the  froe  deceleration 
of  a circular  plate,  obeying  the  Tresca  yield  condition,  are  of  course  well  )cno*m, 
being  due  to  Hopkins  and  Prager  [9,  10]  and  Wang  and  Hopkins  [11].  They  are 
here  cited  as  special  cases  of  the  general  family  of  mode  form  solutions  of  this 
structure . 

(B)  If  we  instead  adopt  Mises'  yield  condition  and  associated  flow  rule, 
the  calculation  of  mode  shapes  and  accelerations  must  be  carried  out  by  a nu- 
merical method.  We  have  done  this  by  several  iterative  schemes  and  by  a finite 
element  program.  Iterative  schemes  are  easy  to  devise  in  cin  intuitive  way,  and 
usually  work  well  in  some  range.  They  have  an  unfortunate  habit  of  failing  to 
converge  outside  of  certain  ranges  of  the  parameters  involved.  When  they  do 
converge,  they  furnish  essentially  exact  solutions  at  small  computer  cost.  Here 
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i 


we  briefly  f'utlin'e  the  iterative  scheme  which  was  found  to  work  well  in  the 
present  problem  over  a wide  range  of  the  relevant  parameters.  Our  finite 
element  solution  was  much  more  expensive.  It  is  mentioned  here  because  in 
subsequent  work  extending  the  mode  technique  to  large  deflections  it  was  found 
to  have  a somewhat  wider  domain  of  convergence  than  the  iteration  scheme. 

The  iterative  solution  of  the  present  problem,  where  only  small  deflections 
are  considered,  was  that  used  in  the  determination  of  mode  shapes  in  a visco- 
plastic plate  [6].  The  use  of  a homogeneous  viscous  representation  of  visco- 

I i 

plastic  behavior  has  been  outlined  above.  Equations  (16)  and  (17)  respectively 

furnish  strain  rates  and  bending  moments  in  terms  of  material  parameters  e^, 

o'  , n'  . The  last  two  are  written  in  terms  of  the  constants  of  Eq.  (I5a)  as 
0 


o'  = po 

O 0 


n ' = vn 


(40) 


where  o , n are  obtained  with  e from  tests  over  a range  of  strain  rates, 
0 0 

and  u , V are  factors  which  enable  the  homogeneotis  viscous  forms  to  match 
test  data  in  a particular  range  of  strain  rates.  We  have  used  the  following 
for  a conservative  and  accurate  matching  with  strain  rate  test  data  [5,  6,  8]: 


1 + B 


1/n 


jl/n 


u = 


1+6 


1/n 


a/^'“ 


(41a ,b) 


where  0 = e /e 
e o 


e being  the  effective  strain  rate  at  which  the  two 
e 


representations  are  matched.  The  function  Eqs.  (18)  serves  as 

in  a plane  stress  problem,  while  serves  as  6 in  a small  deflection 

plate  problem.  As  previously  mentioned,  rigid-perfectly  plastic  behavior  is 
obtained  by  taking  w = 1 and  v very  large,  corresponding  to  putting  6 +■  0 
in  Eos.  (41).  We  outline  the  iterative  scheme  in  the  form  for  a homogeneous 
viscous  behavior,  so  that  either  a rate  sensitive  plate  with  p = 0 (impulsive 
motion)  or  a rigid-perfectly  plastic  plate  with  arbitrary  p can  be  treated. 
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The  equation  of  motion  (28)  may  be  integrated  again  and  put  in  the  form 


■r 

■r 

t 

= 

dr' 

ft,  ^ ♦fdf 

o 

o 

+ (42a) 


where 


P = M ’ - HR— 


(42b) 


and  is  a constant  of  integration.  Now,  from  the  constitutive  equations 

(17),  with  the  mode  solution  form  Eq.  (19)  and  the  curvature  rate  expressions 
Eqs.  (20)«  the  dimensionless  moments  are 


= - 5-  W 


,1/n' 


X R 
0 


(43a) 


O 


where 


X = - [(♦")^  + ^'4"  t (^')^T^^ 
/3  L r r J 


(43b) 


(43c) 


using  Eq.  (43a)  in  Eq.  (42a),  we  write 

r 


A"  — A ' s 

’ 2r* 


|,(mg-mr)dr'-  ipr^  ♦ w^ 


- 3 " 


(44a) 


r R 
0 


l/n'-l 


with 


m.-m 
0 r 


'h[- 

V r 


1/n' 


l/n'-l  I.,,  ^ L^, 

5^  Jr’ 


(44b) 
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We  may  consider  a rigid-perfectly  plate  by  putting  w = 1 and  taking  n'  = vn 
very  large  so  that  the  velocity  disappears  from  the  right  hand  side.  Al- 

ternatively, to  consider  a viscoplastic  material  we  must  put  p = 0 ; again 

* **  * ^ • 

is  eliminated  since  w^/w^  can  be  taken  constant  and  obtained  from 

Eq.  (26).  In  either  case. 


<1> 

= RO?,Cj^) 


(45) 


where  the  right-hand  side  is  known  numerically,  apart  from  the  constant  C,  , if 

<1> 

a first  approximation  ^(r)  is  inserted  in  the  right-hand  sides  of  Eqs.  (41) 

and  in  Eq.  (26).  Integrating  Eq.  (45), 

■r 

r',C^)dr'  (46) 

0 

(An  integration  constant  has  been  taken  as  zero  for  a plate  that  is  continuous 
at  r = 0.)  The  constant  must  be  determined  from  the  edge  condition  at 

r = 1:  for  a simply  supported  plate  the  left-hand  side  of  Eq.  (44a)  vanishes, 
while  for  a clamped  plate  that  of  Eq.  (46)  is  zero.  In  either  case  we  obtain 


<1> 

4'(r)  = S(r)  (47) 

<1> 

where  S is  the  right-hand  side  of  Eq.  (46)  with  the  constant  evaduated 


and 


♦(r)  =■ 


>h> 

S(r' )dr* 


(48) 


-'r 

Finally,  a second  approximation  is  obtained  by  multiplying  the  right-hand 

side  of  Eq.  (48)  by  a factor  which  renders  i(0)  > 1 , giving 


<2> 

♦ (r) 


l(r. 


)dr 


1 

<1> 

S(r')dr' 


r 


(49) 
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<1> 


This  is  put  in  place  of  ^(.r)  and  the  cycle  is  repeated.  Convergence  is  shovm 

* — “ 
by  the  approach  of  for  chosen  p (or  of  w/w^  for  p = 0)  to  a con- 

n 


stant  value,  and  of  the  normalizing  factor 


S(r' )dr'  to  unity. 


' 0 

The  midpoint  acceleration  as  function  of  load  pressure  is  shown  nondimens- 

ionally  in  Fig.  5 for  both  yield  conditions  and  both  types  of  edge  constraint. 

The  shapes  of  the  mode  form  responses  for  the  plates  governed  by  Mises'  yield 

condition  are  illustrated  in  Figs.  6 and  7.  These  evidently  are  close  to  the 

2 

shapes  for  the  Tresca  material,  for  pR  greater  than  about  20. 


Conclusions 

Mode  form  solutions  of  structures  undergoing  dynamic  plastic  deformation 
have  utility  in  approximation  techniques;  they  have  basic  significance 
as  "natural"  response  patterns  to  which  structures  with  arbitrary  starting 
conditions  tend  to  converge.  In  this  paper  we  have  illustrated  sufficient 
conditions  of  loading,  material  behavior,  etc.  for  such  solutions  to  exist. 

We  have  discussed  examples  of  their  calculation  and  have  illustrated  the  links 
between  response  of  dynamically  loaded  structures  and  static  limit  load-plastic 
collapse  analysis.  It  is  hoped  that  the  discussion  will  promote  understanding 
and  use  of  these  properties  of  dynamic  plastic  response. 


h' 
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